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We study the thermal emission of photons from hot and dense strongly interacting hadronic 
matter at temperatures close to the expected phase transition to the Quark-Gluon Plasma (QGP). 
Earlier calculations of photon radiation from ensembles of interacting mesons are re-examined with 
additional constraints, including new production channels as well as an assessment of hadronic form 
factor effects. Whereas strangeness-induced photon yields turn out to be moderate, the hitherto not 
considered i-channel exchange of oj-mesons is found to contribute appreciably for photon energies 
above ~ 1.5 GeV. The role of baryonic effects is assessed using existing many-body calculations of 
, lepton pair production. We argue that our combined results constitute a rather realistic emission 

• rate, appropriate for applications in relativistic heavy-ion collisions. Supplemented with recent 

' evaluations of QGP emission, and an estimate for primordial (hard) production, we compute photon 

, spectra at SPS, RHIC and LHC energies. 
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I. INTRODUCTION 



Electromagnetic radiation (real and virtual photons) from central high-energy collisions of heavy nuclei {A-A) has 
the potential to directly probe the high temperature and high density phases of these reactions. This is so because 
, photons and dileptons suffer little final state interaction, and also because their emission rate is a rapidly varying 
^ function of local intensive variables. It is therefore hoped that electromagnetic (e.m.) observables reveal novel features 
, of strongly interacting systems under extreme conditions, including the quark-gluon plasma (QGP) 0, which is but 
OO ■ one facet of a rich many-body problem. 

I More than a decade of careful experiments and of their theoretical interpretation has made clear that a proof of the 
existence of the QGP requires analyses of many complementary measurements. This, in particular, applies to thermal 
production of photons and dileptons, since their emissivity is encoded in the same theoretical quantity, i.e. the e.m. 
. current-current correlation function in strong- interaction matter |^]. 

In this article we focus on the study of real photon emission (see Ref. Q for recent reviews). An identification 
of the QGP component in the spectra requires a reliable assessment of competing sources, most notably primordial 
production in hard nucleon-nucleon (TV- TV) collisions 0|, as well as the later hadronic emission. Owing to the vanishing 
Clni invariant mass of real photons, all three sources generate essentially smooth spectra, measured as a function of 
^ ■ transverse momentum, qt. The key properties are then their strength together with a typical spectral slope. For 
thermal radiation the latter is governed by temperature (as well as collective matter expansion), whereas the overall 
normalisation is set by a combination of microscopic production mechanisms and the space-time volume of the 
interacting matter. Naively, one expects QGP radiation (from large temperatures and small volumes) to dominate at 
large qt, and hadron gas (HG) emission at small qt. In practice, i.e. within a heavy-ion environment, both ends of the 
• thermal spectra are masked: there is a large background at low qt, and the primordial component from hard N-N 
collisions feeds the high qt. The primordial yield in A- A collisions is expected to be modified (as compared to the 
free N-N case) by nuclear transverse-momentum (fcy-) broadening 0,0 and shadowing 7]. In general, the strength 
of thermal emissivities in both hadron gas and QGP can also be substantially affected by microscopic in-medium 
effects. Furthermore, collective (transverse) expansion induces a spectral hardening being more pronounced in the 
later (hadronic) phases. 

Most of the above issues have been addressed in various forms before, e.g. within hydrodynamic simulations 
transport /cascade models 0, 0, 0| or simple fireball models 0, 0| (where consistency with dilepton observables 
has been emphasised). In the present paper we are mostly concerned with hadron gas emission, improving previous 
analyses in several respects. Trpai-meson gas rates are revisited and extended, including strangeness reactions, heavy 
resonances as well as f-channel exchange of the Ci;(782). The effects of baryons - which are known to be of prime 
importance in the dilepton context 0| (as recently further supported by measurements at lower SPS energy ~ 
extracted based on the same spectral densities that successfully describe the dilepton data, thus ensuring consistency. 
This also holds for our modelling of the space-time evolution of heavy-ion reactions, which, in particular, accounts 
for effects of chemical off-equilibrium in the hadronic phase. Finally, for comparison with data from SPS (and for 
predictions at RHIC and LHC), hadronic emission is supplemented with an estimate of hard (primordial) production, 
an assessment of the Cronin effect, as well as QGP radiation (complete to leading order) within a schematic evolution 
of the collision dynamics. 

The article is organised as follows. In Sec. ^ we briefly recall the definition of emission rates, followed by the 
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evaluation of hadronic photon radiation from (i) a hot meson gas using a massive Yang-Mills (MYM) approach, 
extended to the SU{S) (strangeness) sector, with a quantitative assessment of hadronic form factors; (ii) a baryon- 
rich environment, as well as, (iii) additional meson sources beyond the MYM framework. The total hadronic rate is 
then compared with a recent leading-order QGP result. In Sec. IIIII we employ a simple dynamical model to compare 
integrated yields with experimental measurements at SPS, and to provide predictions for RHIC and LHC. Finally, 
Sec. IIVI contains a summary and conclusions. 



II. THERMAL PHOTON RATES FROM HADRONIC MATTER 



A. Definitions 



In this work we will adopt two different approaches to calculate thermal photon production from hot and dense 
matter, taking advantage of the respective strengths of each one. 

Within the thermal field theory framework, the differential photon emission rate from an equilibrated system can 
be written in compact form as |3] 

^ P,u W^'^ (1) 



where 



P,.=^^a I ^0-^el{\,p) eA\,p) 5''^\p- q) (2) 

denotes the (polarisation-summed) photon tensor to lowest order in the e.m. coupling a. The key quantity is the 
hadronic tensor W^^'^{q] fJ.b,T), which contains the effects of strong interactions to all orders and is directly related to 
the imaginary part of the in- medium photon self-energy (or e.m. current-current correlation function), flcn, via 

IV^" = {-2) f''{qo;T)lmUZ (3) 

with f^{qo; T) = l/(exp[(7o/r] — 1) the Bose distribution function (note the similarity of Eq. to the dilepton rate, 
which is obtained upon replacing with the lepton tensor L^j_^). The photon rate can therefore be cast into the 
form 

- -4 f'^iq^-.T) imnL(go = <z;T) . (4) 

Invoking the vector dominance model (VDM) flem can be directly related to in- medium vector- meson spectral 
functions (see below), which makes Eq. Q particularly suitable for non-perturbative model calculations at low and 
moderate energies and momenta. 

Alternatively, the photon emission can also be computed from relativistic kinetic theory. For a given reaction 
1 + 2^3-1-7, the pertinent rate becomes 



qo- 



dR^ f d^pi d^P2 d^P3 ,4x(4), , , M,„2 /(i?l)/(i^2)[l±/(i?3)] 



d^ - J 2{2nYE, 2{2.fE, 2{2.fE, ^^^^ ^ + P2 ^ P3 + q) \M\ (5) 

where M is the invariant scattering matrix element. For photon self-energies of no more than two-loop order, its 
imaginary part reduces to tree level diagrams, in which case it is usually more convenient to find thermal rates with 
the kinetic formula, Eq. ||SJ). The perturbative character of this approach thus makes it advantageous over Eq. at 
moderate and high energies and momenta. 



B. Meson Gas and Hadronic Form Factors 



To describe photon-producing reactions in a gas consisting of light pseudo-scalar, vector and axial vector mesons 
(tt, K, p, K*,ai) we employ the massive Yang-Mills (MYM) approach, which is capable of yielding adequate hadronic 
phenomenology at tree level with a rather limited set of adjustable parameters. Vector and axial vector fields are 
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implemented into an effective nonlinear cr-model Lagrangian as massive gauge fields of the chiral U(3)lX U(3)fl 
symmetry 0, : 



£ = iF^TiDf.UDt'W + ^F^TtM{U + W - 2) 

-iTr {F^.F^'"' + F^^F"^""") + mgXr (A^A^'^" + A^A^'') + ^iTtF^^U F^^^^W 
-i^Tr {D^UD^^F^^"" + D^U'' D^F^"-") . 



In the above, 



(6) 



U ■ 



exp 



2i (l>iX 



F^ 



V2 



Af:^^{V^ + A^) , 



M = 



L,R 



exp 



2i 



^.A^^'^-^go[A^,^,Ai'^] 



igoA^^U 



igaUAi 



m\ + 



7!^ 



(7) 



Note that F^ — 135 MeV and that is a Gell-Mann matrix. 0, and are, respectively, the pseudo-scalar, vector 
and axial vector meson matrices. Note that this form of the interaction enables a coherent treatment of the strange 
and nonstrange fields. 

Let us first concentrate on the nonstrange sector (Trpai). Following Ref. 20], the four unspecified parameters in 
the Lagrangian can be inferred from the masses and widths of the p and ai, allowing for two possible solutions: 



(I) : g = 10.3063, 7 = 0.3405, ^ = 0.4473, mo = 0.6253 GeV ; 

(II) : g ^ 6.4483, 7 = -0.2913, C = 0.0585, too = 0.875 GeV 



(8) 



where g = 5o/\/l ~ 7- I'^ the absence of additional empirical constraints, the use of one set over another is difficult to 
justify. However, in Ref. |2H it has been suggested to invoke the experimental determination of D- to S'-wave content 
in the final state of the Oi — > pir decay. For the two parameter sets given above one finds 



(I) : D/S = 0.36 

(II) : D/S = -0.099 , 



(9) 



to be compared to the experimental value |22l | — 0.107± 0.016. From here on, we therefore employ parameter set II. 

We proceed with a systematic evaluation of all processes generating photons based on the interaction vertices 
contained in Eq. The explicit reactions considered include all possible s-, t- and u-channel (Born-) graphs for 
the reactions: X + Y^Z-\-^, p^Y-\-Z-\-^ and K* Y + Z + 'j. For X, Y, Z we have each combination 
of p,TT, K* , K mesons which respect the conservation of charge, isospin, strangeness and G-parity defined for non- 
strange mesons. The axial ai meson has been considered as exchange particle only (the ai — > decay is automatically 
incorporated via s-channel irp scattering). Using Eq. (|SJ| the thermal photon production rates are readily obtained 
from the coherently summed matrix elements in each channel, and convenient parametrisations thereof are given in 
Appendix ^ The parametrisations for the nonstrange reaction channels in the present work differ from the ones 
quoted in Ref. J3j, which are based on Ref. 20|, in two respects. First, in the latter article the amplitude in two 
channels, as written, violates the Ward identity |2^. Second, the choice of parameters underlying Ref. yields an 
D/S ratio which is at variance with the experimental value. This is corrected here, together with amplitudes which 
have been verified for current conservation. As in Ref. [23|, parameter set (II) yields smaller emission rates than set 
(I) for the (leading) irp ttj process at high energies. In addition, our results using set (II) are another ^40% smaller 
than the corresponding ones in Ref. . 

An important element in applying effective hadronic models at moderate and high momentum transfers is the use 
of hadronic vertex form factors simulating finite-size effects, which are not accounted for in the Trpoi calculations 
of Ref. [23 ■ For some of these reactions form factor effects have first been studied in Ref. and found to give a 
typical net suppression over the bare graphs by an appreciable factor of ~3 at photon energies around go— 2-5 GeV. 
We perform an estimate of their effect in the present context as follows. To be consistent with the procedure adopted 
before for dilepton production ,2^, we assume a standard dipole form for each hadronic vertex appearing in the 
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FIG. 1: (Color online) The effect of hadronic form factors on a typical photon-producing reaction in a nonstrange meson gas, 
TTp 7r7, for 3 different temperatures. Solid and dashed lines are without and with the inclusion of F{t)*', respectively. 



amplitudes, 

We then approximate the four- momentum transfer in a given i-channel exchange of meson X by its average t according 
to 

\2K^-t) ~iE^ Jo {m\-tY{2K-^ -tf ' ^^^^ 

The averaging procedure allows to factorize form factors and amphtudes which much facilitates the task of rendering 
the final expression gauge invariant, since we can then simply multiply the rate parametrisations with F(t)^ employing 
a typical hadronic scale, A = 1 GeV |2^, for the cutoff parameter. The impact of the form factors is indicated by 
the difference of solid and dashed curves in Fig. ^ which should still be considered as a conservative estimate of the 
suppression effect. The reduction of the rate in the 2-3 GeV region amounts to a factor of ^i, quite in line with the 
earlier estimates of Ref. '25^ . 

Rather little attention has been paid in the literature so far to the calculation of photon emission rates involving 
strange particles, mostly because existing analyses have found them to be quantitatively suppressed. In Ref. |23|, the 
channel Ki was investigated and shown to be appreciable relative to non-strange sources only in a limited 

kinematical domain. Also, results obtained with SU{3) chiral reduction formulae p^ . coupled with an expansion in 
temperature, suggest that the strangeness contribution is not large. Here, we seek to quantify the latter relative to 
the irpai emissivities within the same effective Lagrangian framework as encoded in the SU(3) extension implicit in 
Eq. To optimally reproduce the (measured) hadronic phenomenology, we are, however, lead to decouple the non- 
strange axial vector meson (ai) from the strangeness sector. This allows to simultaneously satisfy the (electromagnetic) 
Ward identities and fix both the strange vector mass, mK*=8Q5 MeV, independent of the p mass, and the universal 
coupling constant as to match the empirical value 2^ of the K* width, T{K* — > _ftr7r)~50 MeV. The specific channels 
included are: ttK* — > K-y, Kp Ky, ttK K*'-f, KK* nj, KK — > p7, and K* ttKj. Again, all possible 
isospin combinations are accounted for in the rate calculations, which have been parameterized in functional form in 
Appendix^ For a temperature of 200 MeV, the leading production channels are shown in Fig.|21 with hadronic form 
factors implemented following the same procedure, Ea. (|ll|) . as before (the last two contributions enumerated above 
have been omitted, as they represent negligible increments). At all energies of practical relevance (go > 0.5 GeV) and 
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FIG. 2: (Color online) Photon-producing reaction rates involving strange mesons at a temperature r=200 MeV with form 
factor effects included. 



including form factors, the K*ii —f reaction, mediated by t-channel 7r-exchange, turns out to be the main emission 
source, which is in complete analogy to the 7rp — > 7r7 reaction in the nonstrange sector. In line with estimates in (29| . 
the former constitutes ~40% of the latter around qo=l GeV, being reduced to about 20% at go=3 GeV. 



C. In-Medium Selfenergies with Baryons 

It is important to realize that thermal emission rates of dileptons and photons are intimately connected, both 
being based on the e.m. current-current correlator, albeit evaluated in distinct kinematical domains, i.e. timelike 
(M2 = ql-q^ :>{)) vs. HghtHke (M^ = 0), respectively. The latter may imply the prevalence of different processes in 
the corresponding observable, but consistency can and should be tested. In particular, we recall that baryonic sources 
are very important 0, [s^l for understanding the observed excess in low-mass [M < 1 GeV) dilepton production in 
(semi-) central Pb-Au collisions at both full (160 AGeV) [3l| and lower (40 AGeV) SPS energy. It is therefore 
mandatory to scrutinize the role of baryons in photon production, especially since most investigations thus far not 
revealed substantial contributions [s^, IH, 0| • 

We here make use of the hadronic many-body calculations of the in- medium p(770) spectral function [ssl Is^ Is^. 
which, when evaluated for 0, directly yield pertinent photon emission rates via Eq. Within the vector 

dominance model (VDM), one has (schematically) 

4 

ImHen, - ^Im^V (12) 

(my, gv and ImDy: vector- meson masses, coupling constants and spectral functions, respectively). In the following we 
focus on contributions arising from the p-meson, which are dominant since g^p/ g^^lQ. In addition to in-medium effects 
in the pion cloud of the p-meson (encoded in a modified two-pion decay width, YipT^j^)^ resonant p-h {h=Tr, K, p, TV, A,...) 
interactions are incorporated through self-energy expressions of type 

S;,^(90,9;T) = I ^-L_[f\u;,{p))-r\u;,{p) + qo)] M^pl{p,q) , (13) 

where the isospin averaged p scattering amplitude Aiph is integrated over the thermal distribution f^{LOh{p)) — 
[exp(w;i(p))/T ± 1]^^ of the corresponding hadron species h with LOh{p) — ymf^+fp. The advantage of writing the 
self-energy in terms of the forward scattering amplitude is that in- medium resonance widths, accounting for higher 
order effects in temperature and density, are readily implemented without facing problems of double-counting. The 
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latter becomes more difficult to keep track of when evaluating higher order topologies in the kinetic-theory a ppr oach 
represented by Eq. ||SJ) . All of the resonances used in constructing the p self-energy are enumerated in Refs. |26ll36l| . 
which also contains more details on how the interactions are constrained by hadronic phenomenology. 

The results from the hadronic many-body approach are compiled in Fig. |21 for two temperature-density values 
characteristic for meson-to-baryon ratios at full CERN-SPS energy (160 AGeV). The solid curve is the net photon 
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FIG. 3: (Color online) Thermal photon production rate (under conditions resembling CERN-SPS P6(158AGeV)+P& collisions) 
in the hadronic many-body approach of Refs. |2^l35ll3^ based on an in-medium p spectral function. Left panel: for temperature 
and baryon chemical potential {^iB,T) = (340, 150) MeV, right panel: {^iB,T) = (220,200) MeV. 

spectrum obtained by taking the full (p- meson) spectral density to the photon point, whereas the long-dashed curve 
represents the non-baryonic (sub-) component. One observes that the low-energy regime, <1 GeV, of the photon 
emissivity is dominated by baryonic effects (quite reminiscent to what has been found for low-mass dileptons) . These 
are mostly due to direct pN resonances such as A(1232), A^(1520), as well as A(1232)A^~^ and NN~^ excitations in 
the two-pion cloud of the p (which, to leading order in density, correspond to t-channel one-pion exchange (OPE) in 
processes of type ttN 'jN) . These contributions should be rather reliable for baryon densities up to at least normal 
nuclear matter density, po = 0.16 fm~^, being constrained by photo-absorption spectra on nucleons and nuclei [3^ 
(including hadronic vertex form factors with rather soft cutoff parameters around 600 MeV). At comparable baryonic 
densities this approach yields about a factor of two more photons than the results obtained in Ref. [s^, where only 
nucleonic degrees of freedom where accounted for within the (on-shell) chiral reduction formalism (see Ref. [39| for 
an update including A(1232) and A^(1520) resonances). 

Beyond ~ 1 GeV, mesonic (resonance) states become the dominant source of photons in the many-body approach, 
which includes radiative decays of w(782), /ii(1170), ai(1260), /i(1285), 7r(1300), a2(1320), w(1420), w(1650), K*{892) 
and ifi(1270). In particular, the cj — + 7r7 decay exhibits an appreciable low-energy strength, consistent with the early 
results of Ref. 25] . Note that all hadronic vertices carry (dipole) form factors with typical cutoff parameters of around 
1 GeV, as extracted from an optimal fit to measured hadronic and radiative branching ratios within VDM [26j : t- 
channel exchange processes between mesons as discussed in Sec. Ill Bl (e.g. OPE or ai-exchange in np tt'j) are not 
implicit in the spectral densities leading to the the results of Fig. They are mostly relevant at photon energies 
beyond 1 GeV and therefore do not significantly figure into bulk (low-mass) dilcpton production, the latter being 
dominated by (transverse) momenta qt < M. As mentioned above, the underlying VDM coupling to the photon 
exclusively proceeds through the p(770) which implies that the strength in the pertinent e.m. correlation function 
beyond mass/energy scales of ~ IGeV is no longer correctly saturated, as it is restricted to two-pion-type states. The 
construction of the total emission rate will be discussed in the following Section. 



D. io i-Channel Exchange and Total Rate 



Before combining the various contributions to the thermal photon rate we investigate one more process of potential 
importance which is not present in the above and, to our knowledge, has not been addressed before. Within the 
SU{2) flavour symmetry, the u; is a chiral singlet, but is known to exhibit a large coupling to 7rp and thus, via VDM, 
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FIG. 4: (Color online) Compilation of thermal photon production rates from hot and dense hadronic matter computed in 
the present work at temperature r=200 MeV and baryonic chemical potential /is =220 MeV (translating into total-pion to 
net-baryon ratio of -^5). Dashed and dotted lines correspond to the non-/strange MYM meson gas emissivities of Sec. Ill B1 
using the parametrisations given in Appendix^ solid line to the p spectral function approach including baryons, and the 
dashed-dotted line is solely due to uj t-channel exchange in vrp ^ 717. 



to 7r7 states. Its s-channel decays have indeed been calculated as early as in Ref. 25] (and are included above), 

but its t-channel exchange in the reaction -np — > 7r7 has not. We here have calculated the pertinent contribution to 
the thermal emission rate using Eqs. jSJ) and with the same coupling and form factor type as for the s-channel 
graph [2^ (corresponding to Fig.|3), see below. 

In Fig. 0] we summarise our results for the thermal photon emissivities from hadronic matter as evaluated in the 
preceding Sections. At low energies, go ^ IGeV, the emission rate from the hadronic many-body approach (p spectral 
function) j36| . with major contributions from baryonic sources, dominates. Between energies of 1 and 2 GeV, meson 
gas emissivities become competitive and eventually dominate the rates at high energies. Remarkably, the w t-channel 
exchange in tt/? 717 is the single most important process beyond energies of go— 2 GeV. The strangeness component 
in the production rate does not exceed 10-15% at any energy. 

We finally have to address the question of how to combine the various hadronic sources, computed in two different 
frameworks (cf. Sec. Ill A|) . into the total emission rate. Two issues arise when simply adding all of the emission rates 
shown in Fig.^] double-counting and coherence. The ai s-channel graph is present in both p spectral function and the 
MYM framework. We remove it from the former, where it plays a minor role, whereas it induces significant interference 
effects in the Ttpoi complex. If coherence is unimportant, t-channel contributions can be evaluated separately. It was 
verified that this was the case for the uj exchange, so that the incoherent addition of the ^-channel contribution is 
justified. 

We believe that it is fair to say that the enumeration of hadronic photons sources given in this Section, together 
with form factor inclusions, currently represents the most realistic evaluation of the full hadron gas emissivity. 



Before turning to applications in heavy-ion reactions, our estimates for hadronic production rates are confronted 
with the ones from QGP emission, in particular with the simple lowest-order HTL-corrected pQCD result |25i||. 



and (a parametrization of) the complete leading-order (in aag) analysis ||4^J, cf. Fig.|Sl Clearly, due to the approx- 
imations implied by each curve, none of them can be expected to be accurate under conditions arising in practice, 
i.e. in the phase transition region. Nonetheless, the observation that the complete leading-order QGP calculation 
(dashed-dotted curve) is similar to the full hadronic result (sum of solid and long-dashed curves) within a factor of 




E. Comparison to QGP Emission 




(14) 
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~ 2 over essentially all (relevant) energies below 3 GeV might not be a mere coincidence. A similar behaviour has 
been found before for dilepton production rates iSQ], perhaps suggesting a type of "quark- hadron duality" for e.m. 
emission close to the expected phase boundary. It would be most valuable to shed further light on this issue from 
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FIG. 5: (Color online) Comparison of HG and QGP photon production rates at r=200 MeV. Solid line: hadronic many-body 
approach of Refs. l2 a . [Sa. l3ql (solid curve), dashed line: mesonic contribution including hadronic form factors, dotted line: 
simple pQCD result j25(| according to Eq. 11411 . dashed-dotted line: complete leading-order QGP emission |40| . 

first principle lattice calculations, which, at the moment, are only reliable at sufficiently large invariant masses where 
they are, as expected, close to perturbative results [4l| . 

III. PHOTON SPECTRA IN ULTRARELATIVISTIC HEAVY-ION COLLISIONS 

A. Hard Photons and Thermal Fireball Evolution 

For a realistic comparison with direct photon spectra as extracted in heavy-ion collisions two further ingredients 
are required. 

First, the thermal rates of the previous sections have to be convoluted over the space-time history of the reaction. 
Assuming that thermal equilibrium can be established and maintained, hydrodynamic simulations are the method of 
choice, see e.g. Ref. 10]. Here we employ a more simple fireball model |3a.l4a |. which incorporates essential elements 
of hydrodynamic calculations. The fireball evolution is started at a "formation" (or thermalization) time tq < 1 fm/c, 
which relates to the initial longitudinal extent of the firecylinder as Az ~ Aj/tq with Ay ~ 1.8 corresponding to the 
approximate rapidity coverage of a thermal distribution. The subsequent volume expansion, Vp b (t) , is carried through 
QGP, mixed and hadronic phases until "thermal" freezeout at r/o=100-120 MeV, where hadrons cease to interact. 
The equations of state (EoS) for QGP and HG are modelled via thermal quasiparticles and a resonance gas (including 
about 50 species), respectively. Based on the conservation of net baryon-number, Nb, and total entropy, S, one is able 
to extract the temperature and baryon chemical potential at any given (proper) time, thereby defining a trajectory in 
the hb-T plane. The transition from the QGP to HG phase is placed at "chemical freezeout" points extracted from 
hadron ratios in experiment |43|| . Consequently, in the HG evolution from chemical to thermal freezeout, hadrons 
stable under strong interactions (pions, kaons, etc.) have to be conserved explicitly by introducing associate chemical 
potentials (/x.^-, jiK, etc.). This has not been done in previous calculations of thermal photon production 8, 91. llOlj . 
and induces a significantly faster cooling in the hadronic phases ,29j. In addition, at collider energies (RHIC and 
LHC), the conservation of the observed antibaryon-to-baryon ratio (which at midrapidities is no longer small) in the 
hadronic evolution becomes important . An accordingly introduced (effective) chemical potential for antibaryons 
has been shown to impact the chemistry at later stages appreciably |44| (in particular, it is at the origin of large 
meson-chemical potentials, again implying faster cooling). For each collision energy, the value of the specific entropy, 
S/Nb, is fixed to reproduce observed hadron abundances. The total yield of thermal photons in an A- A collision then 
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follows as 



dN. 



th 



ermal 



g° '"rf3g fa) = ^ I JdrVFBir) (^90^) , (15) 

averaged over a rapidity interval, [ymin,ymax] according to the experimental coverage {Ay — ymax — ymin)- To 
incorporate transverse expansion in the spectra, in the thermal rest frame isotropic photon momentum distributions 
are boosted into the laboratory (lab) frame using an average of about 70% of the time-dependent transverse (surface) 
expansion velocity at each moment in the fireball evolution. 

Second, an additional contribution to direct photon spectra arises from prompt photons in primordial N-N collisions. 
The minimal baseline for a heavy-ion reaction constitutes the collision-number scaled expectation from proton-nuclcon 
collisions. An accurate description thereof is still a matter of debate 0|, so that we here employ empirical scaling 
relations in xt = 2qt/^/s, extracted from fits to data in Ref. For xt > 0.1, corresponding to fixed target energies 
{■\/s < 50 GeV and photon transverse momenta qt above 2 GeV), the cross section fit (at midrapidity, y—Q) reads 

d'aPP _ i^sr-^ pb 
'^"^ - feF^ ' 

whereas for xt<0.1, corresponding to collider energies {y/s > 200 GeV and <Zt<10 GeV), 

The naive extrapolation to a collision of two nuclei A and B at impact parameter b predicts the prompt photon 
spectrum to be 

(l]\f prompt (PaPP 

90^5 {b;quy = Q;V^) = qo-j^ ABTAsib) 

dPq d^q 

= q^^ — 18 

d^q a- ^ ' 

with T4B: nuclear overlap function, Ncoii- number of primordial N-N collisions, cr™: inelastic N-N cross section (we 
also supplement Eq. ()18|l with a smooth cutoff for go < 2GeV, where the parametrization (|16|l is no longer reliable). 
The lack of a consistent microscopic description of photon production in p-p complicates the task to assess nuclear 
corrections, such as shadowing or intrinsic kx broadening (Cronin effect), see e.g. Ref. |3] for a recent discussion. As 
a substitute for a more rigorous calculation, we here adopt the following strategy: since the intrinsic kx effects at the 
N-N level are in principle contained in the parametrization, Ea. (|16|l . the nuclear effect is approximated by fitting an 
additional (nuclear) fcy-smearing to p-A data. The latter is modelled by folding the parameterized spectrum over a 
Gaussian distribution 

m) = e-^/^-^) . (19) 

The result, together with proton-nucleus data on photon production is shown in Fig.El The data have been scaled 
to the 10% central Pb-Pb cross section at 158 AGeV according to the procedure used in Ref. Fitting the 

p-A single photon data in this fashion, an adequate reproduction of the experimental measurements emerges with 
(Afc2)~0.1-0.2 GeV^. 

In principle, there is a third source of photons corresponding to emission after initial nuclear impact, but before the 
formation time tq (the " pre-equilibrium" contribution). It is difhcult to assess both theoretically and experimentally; 
a rough (but uncontrolled) estimate might be had by choosing a somewhat smaller formation time. Note that, in 
principle, the modelling of those contributions is accessible to ab-initio simulations [49l |. 



B. SPS 



In this section we compute transverse momentum spectra at midrapidities from 10% central P6(158AGeV)-|-P6 
collisions for which photon spectra have been measured by WA98 ,4^]. Let us first focus on thermal emission; QGP 
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radiation is always calculated with the complete leading-order result hadronic radiation as the sum of the 

hadronic many-body H^H^, TrpK*K gas (within MYM) and uj i-channel exchange contributions with form factors 
as discussed in Sec. Ill Dl 

In Fig. we display results for a rather standard fireball evolution with initial zo=1.8 fm (corresponding to tq ~ 
1 fm/c and initial temperature Ti — 205 MeV) and final temperature T/o— 110 MeV reached after a total lifetime 
of ^ 13 fm/c. Due to transverse expansion, the total hadron gas yield outshines QGP emission at all momenta. 
This is in close reminiscence to the calculations of intermediate dilepton spectra within the same framework |42j |. 
where QGP radiation was found to constitute about 30% of the thermal component that was able to reproduce the 
excess observed by the NA50 collaboration js^ (see also Ref. ^1]). As expected, photons of baryonic origin prevail 
in the spectrum for gt<l GeV; this region is thus intimately related to the low-mass (and low-transverse momentum) 
dilepton enhancement observed by CERES/NA45 [sj]- The same feature has been found 0] when comparing the 
hadronic many-body contributions to the upper limits in S{200AGeV)+ Au by WA80 [s^- In the right panel of Fig.[7| 
we illustrate the sensitivity of the hadronic thermal emission to properties of the fireball evolution. When increasing 
the thermal freezeout time from our (SPS) default value of 106 MeV to 135 MeV, the yield at gt<lGeV is reduced 
by up to 30%, whereas it is essentially unchanged beyond qt—2 GeV, thus reflecting emission close to Tc. On the 
contrary, if the additional boost on the photons due to the transverse expansion is neglected, the high-momentum 
spectrum is reduced appreciably (by a factor of ~3 already at qt=2 GeV), whereas the low-momentum region is only 
mildly affected. 

As is well-known, high-energy (-mass) photon (dilepton) emission is rather sensitive to initial temperatures in 
heavy-ion reactions, due to the large (negative) exponents in the thermal factors. This is confirmed by our results 
for the QGP contribution in Fig. |S1 when decreasing the formation time from 1 to 0.56 fm/c, the latter implying 
Ti— 250 MeV. Another effect that has been ignored in available hydrodynamic calculations so far is associated with 
corrections to the QGP equation of state. The standard assumption is that of an ideal (massless) gas with an effective 
number of flavours Nf = 2.5 corresponding to a total degeneracy dQGP=(10.5A^/-|-16)=42 (in our default calculations 
we use dQGP=40). However, lattice gauge theory results [S^l indicate ~20% smaller values than the ideal gas for the 
thermodynamic state variables in the for SPS energies relevant temperature region r=l-2 Tc- Implementing such a 
reduction into the entropy density (which is the relevant quantity for the fireball evolution) by using donp—32 yields 
an increase of the initial temperature from 250 MeV to 270 MeV without decrease in initial volume |29l Is^ l . The 
resulting thermal photon spectrum from QGP radiation shows appreciable sensitivity to this modification at high 
energies, cf. solid vs. dotted curve in Fig. |S1 This sensitivity does not persist into the low-energy region 90 ^ 1 GeV, 
and thus does not affect low-mass dilepton production [33 |. 

Let us now turn to a comparison with the recent measurements of WA98 48:j. Our baseline scenario consists 
of thermal emission (hadronic and QGP) from the expanding fireball with Ti=205 MeV, supplemented by prompt 
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FIG. 6: (Color online) Direct photon data in proton-Carbon collisions, scaled to central Ph-Ph collisions at SPS energies (see 
text for details). The curves show the effect of the broadening of the primordial photon spectrum generated by the nuclear 
medium. The data are from Refs. ^ (E629) and [i^] (NA3). 
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FIG. 7: (Color online) Integrated photon emission from various thermal sources in an expanding fireball model for central 
Pb+Pb collisions at SPS. Left panel: hadronic emission (T < Tc=175 MeV) from the meson gas component (dashed line) and 
the in-medium p spectral function (solid line) compared to QGP emission {Tc < T < Ti— 205 MeV) (dashed-dotted line). Right 
panel: Sensitivity of the total HG yield to thermal freezeout (long-dashed line: Tfo~lS fm/c corresponding to r/o=106 MeV, 
dotted line: r/o=9 fm/c corresponding to Tfo~lS5 MeV) and transverse flow (dashed-dotted line: Tfo~lS fm/c with no 
transverse boost of the emission source). 
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FIG. 8: (Color online) Sensitivity of the QGP photon emission yield within an expanding fireball model for central Pb+Pb 
collisions at SPS. 



(pQCD) photons from primordial N-N collisions without any nuclear effects, cf. Fig.^ Up to transverse momenta of 
about 1.5 GeV the data (upper limits) are essentially saturated by thermal radiation from the hadronic phase. This 
is gratifying to note since this regime, as discussed above, is directly related to the low-mass dilepton excess observed 
by CERES/NA45 jHlp, which can be successfully described within the same approach [s^. Beyond 3 GeV, prompt 
photons dominate, but do not seem to provide enough yield to account for the data. Since the hadron gas emission 
is essentially fixed and describes well the low-energy regime, three possibilities are left for the origin of discrepancies 
above go— 2 GeV: (i) modifications of the prompt yield, (ii) pre-equilibrium emission, (iii) larger QGP radiation. In 
the following, cases (i) and (iii) (or a combination thereof) will be investigated. 

First, we study the effects of the initial temperature on the photon spectrum. The exercise in the right panel of 
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FIG. 9: (Color online) Thermal plus prompt photon spectra compared to data from WA98 E^l for central Pb+Pb collisions at 
SPS . 

Fig. 13 is repeated adding all sources discussed here, cf. Fig. EB Clearly, when going to (for SPS conditions) rather 
short formation times of tq — 0.5 fm/c, coupled with non-perturbative (suppression) effects in the QGP EoS, a rather 
good reproduction of the entire spectrum can be achieved. This statement agrees with the hydrodynamic analyses of 

Refs. la 113 • 
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FIG. 10: (Color online) Effects of various initial temperatures on the total photon spectra in Pb+Pb collisions at SPS, compared 
to data from WA98 ||4a|. 

The second possibility relates to the nuclear Cronin enhancement, which we implement as outlined in the previous 
Section. The usual assumption to extrapolate nuclear broadening effects on e.g. 7r° or 7 spectra is that 

{Akl)AA = N{Akl)pA , (20) 

with N=2 [53. Alternatively, based on a careful analysis of the target (A) dependence in p-A collisions, it has been 
suggested in Ref. |56j | that the Cronin effect is due to no more than one semi-hard collision prior to the hard scattering. 
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and therefore saturating as a function of the A^'-A'' colhsion number. In this eventuahty, N < 2. Recalhng that, from 
Fig. ini (A/c^)=0. 1-0.2 GeV^ gives a reasonable description of the 7 spectra in p-C, (AA:|,)-values between 0.2 and 

0. 3 GeV^ seem appropriate for central Pb-Pb collisions. One should also note that the pertinent spectral enhancement 
in the qt—S GeV region amounts to a factor of around 3, which is quite consistent with the nuclear enhancement in tt'^ 
production observed in the same experiment js^. In Fig. II II we have combined the baseline thermal yield {tq=1 fm/c, 

1. e. 7i=2G5 MeV) with 3 values for the nuclear /cT-broadening, i.e. {Ak'^)—0, 0.2 and 0.3 GeV^. The thermal plus 
Cronin-enhanced pQCD spectra provide good description of the WA98 data, even with an initial temperature as 
low as Ti=205 MeV. This constitutes one of the main results of our work: the photon spectrum in nucleus-nucleus 
collisions at SPS energies is perfectly compatible with "moderate" initial temperatures. It also complements, within 
a common thermal framework, earlier descriptions of low- and intermediate- mass dilepton spectra [3^ E3| , as well 
as J/tp an d tp ' production systematics [s^, as observed by the CERN-SPS experiments CERES/NA45 0,1^ and 
NA5G iSlia, respectively. 
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FIG. 11: (Color online) Effects of the nuclear broadening of the primordial photon spectrum on the measured spectrum. All 
the sources discussed in this paper are included in the space-time evolution. 

It is also of interest to quote the values of the transverse momentum where the pQCD yield exceeds the total 
thermal one; these are qt=2.55, 1.7 and 1.55 GeV, corresponding to (Afcf,)=0, 0.2 and 0.3 GeV^, respectively. Again, 
this compares well with the calculation of intermediate-mass dileptons in Ref. [i^ , where the Drell-Yan contribution 
was found to exceed the thermal one at M^^~2 GeV (note that the Cronin effect is expected to be less pronounced 
for quarks than for gluons, and thus will affect the Drell-Yan process less than prompt photons). 



C. RHIC and LHC 



At collider energies the space-time evolution of the expanding QGP and hadronic fireball is expected to change in 
several respects. First, higher energies entail larger charged particle multiplicities per unit rapidity, dNch/dy. In central 
Au+Au collisions at full RHIC energy {^/s = 200 AGeV) experiments have found |60ll6ll| about a factor of 2 increase 
as compared to maximum SPS energy {^/s = 17.3 AGeV). Extrapolations into the LHC regime {^/s — 5500 AGeV) 
suggest another factor of up to ~4 enhancement over the RHIC results. 

Second, the net baryon content at midrapidity decreases, implying small baryon chemical potentials at chemical 
freezeout, e.g. — 25 MeV at RHIC-200. At the same time, the observed production of baryon-antibaryon pairs 
strongly rises, resulting in total rapidity densities for baryons at RHIC that are quite reminiscent of the situation 
at SPS energies '62]. This observation not only necessitates the explicit conservation of antibaryon- number between 
chemical and thermal freezeout j44j (see above), but also requires to evaluate baryonic photon sources with the sum 
of the baryon and antibaryon density (strong and e.m. interactions are CP-invariant). 

Third, the transverse expansion (i.e. flow velocity) increases by about 20% from SPS to RHIC (presumably further 
at LHC), whereas the total fireball lifetime does not appear to change much. The latter, however, is likely to increase 
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at LHC, due the significantly larger system sizes towards thermal freezeout. 

All these features are readily implemented [i^. Ig^ I into the thermal fireball description employed for SPS energies 
above. In addition, the primordial pQCD component changes its xt-scaling behavior which is accounted for by 
replacing the parametrization Eq. (|16|l by Eq. (|17l) . For simplicity we here refrain from introducing a nuclear kT 
broadening, which is expected to be much less pronounced (and/or compensated by shadowing corrections) at collider 
energies. First data on high-pt hadron production in d-Au collisions at ^—200 GeV |3|6^|63 indeed indicate only 
a comparatively small enhancement of around 20-30% over the spectra measured in p-p collisions. 

Our photon predictions for full RHIC energy are summarized in Fig. 1121 The thermal component has been evaluated 
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FIG. 12: (Color online) Integrated photon emission spectra from central Au+Au collisions at RHIC. Short-dashed line: pQCD 
photons from primordial A^-A'' collisions; dashed-dotted line: thermal QGP radiation; long-dashed line: thermal hadron gas 
emission; solid line: total direct photon yield. 

with a typical formation time, tq = 1/3 fm/c, as used before in dilepton and charmonium j58j | applications (it 
is also consistent with hydrodynamic approaches that correctly reproduce the elliptic flow measurements which are 
particularly sensitive to the early phases, see Ref. for a recent review). One notes that the spectrum decomposes 
into essentially 3 regimes: at low energies, qo<l GeV, the major source are still thermal hadrons, whereas at high 
energies, go > 3GeV, prompt pQCD photons dominate. The intermediate region, l<qo<3 GeV, appears to be a 
promising window to be sensitive for thermal QGP radiation. The latter has been calculated assuming chemically 
equilibrated quark- and gluon-densities throughout. It is conceivable, however, that the early QGP phases are gluon- 
dominated, i.e. with quark fugacities much smaller than one (even the gluon fugacities could be reduced). In this 
case, on the one hand, the photon emissivities at given temperature are severely suppressed. On the other hand, if 
most of the total entropy is produced sufhciently early, smaller fugacities imply larger temperatures, thus increasing 
the photon yield. The interplay of these effects has been studied for dilepton production in Ref. 63], where it 
has been found that the net effect consists of a slight hardening of the QGP emission spectrum with a pivot point 
at M~3 GeV. For photons the situation might be even more favorable due to the participation of gluons in their 
production {e.g. g + q ^ gj, as opposed to leading-order qq ee for dileptons). 

We finally turn to the LHC, cf. Fig.^J According to our estimates, assuming a formation time of 0.11 fm/c (trans- 
lating into ri~850 MeV for dNch / dy~3000) , the QGP window extends significantly further in transverse momentum 
than under RHIC conditions, cf. left panel of Fig.^1 although this feature is sensitive to: (i) the formation time (ii) a 
possible chemical undersaturation of the QGP, (iii) nuclear effects on the initial pQCD yield. The transition from HG 
to QGP dominated emission occurs again close to qt=l GeV. In the right panel of Fig.^Jwe illustrate the sensitivity 
of the thermal spectra with respect to the produced charged particle muliplicity within our schematic fireball evolution 
model. For simplicity, we assumed the same formation time and expansion parameters for both Nch—-^000 and 2000. 
We then find that the total integrated photon yield {i.e., for transverse momenta above 50 MeV) scales as A''^ with 
a~1.4 if the same thermal freezeout temperature Tfn^dO MeV is imposed. This value for a is somewhat larger than 
the 1.2 found in hydrodynamic caclulations of Ref. |68j . but confirms the deviation from the naive quadratic behavior 
(a=2). However, the latter is approached (and even exceeded) for the yield at higher transverse momenta; e.g., when 
integrating over qt with a lower bound of 1 GeV (2 GeV), a increases to ~1.9 (2.3). 
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FIG. 13: (Color online) Left panel: integrated photon emission from various thermal sources in central Pb+Pb collisions at the 
LHC; line identification as in Fig. 1131 Right panel: sensitivity of thermal emission spectra from the hadron gas (solid lines) and 
the QGP (dahed lines) to the charged particle rapidity density. For HG emission at Nch='2000, the upper line corresponds to 
assuming the same freezeout time (r/o— 25 fm/c) as for Nch=3000, whereas the lower line corresponds to assuming the same 
freezeout temperature (r/o~90 MeV) (implying a 10% reduced lifetime) . 



IV. SUMMARY AND CONCLUSIONS 



In the present article we have attempted an improved evaluation of hadronic thermal emission rates for real photons, 
suitable for realistic applications in relativistic heavy- ion collisions. In what we think is a more complete treatment 
than has been achieved before, our main findings are: 

(i) Revisited meson gas emissivities built upon an effective Lagrangian of the massive Yang-Mills type lead to about 
40% reduced rates in a irpai gas as compared to previous analyses. An inclusion of strangeness-bearing channels 
has revealed that the latter contribute at the 20% level. A quantitative evaluation of hadronic form factors has 
been performed throughout, mandatory for applications. 

(ii) Photon rates from the baryonic sector have been obtained from a limiting procedure where in-medium p spectral 
densities were carried to the photon point. This procedure makes consistent the real photon analysis with that of 
low invariant mass dileptons, thereby elucidating the role of baryons in photon emission during nuclear collisions. 
Their contributions have been shown to be substantial for photon energies qo<l GeV. 

(iii) As the single most important process at high energies we have identified uj t-channel exchange in the irp ttj 
reaction, which had not been considered before. 

The total hadronic emissivity has been compared to a recent complete leading-order (in strong and electromagnetic 
couplings) QCD calculation for the QGP. In the vicinity of the expected phase boundary both rates turned out to be 
very similar, at all energies of practical relevance. 

The net rates have been folded over a fireball evolution of nuclear collisions. This approach, albeit schematic, 
is consistent with observed hadrochemistry and hydrodynamic expansion characteristics, as well as dilepton and 
charmonium data measured at SPS energies. Using a comprehensive fit of photon cross sections in p-p and p- 
p interactions, an estimate of the Cronin effect (nuclear fcx-broadening) in p-A collisions was first extracted, then 
generalised to central Pb-Pb collisions to address the WA98 photon measurements at the SPS. Combining the complete 
set of hadronic rates with QGP emission and our Cronin-effect estimates on the primordial photon component, we 
are able to reproduce the WA98 data, with moderate values of initial temperature (Ti~200-240 MeV) and transverse 
momentum broadening ((Afc|,)~0.2-0.3 GeV^). Predictions for photon spectra to be measured at higher colliding 
energies (RHIC and LHC) have also been made, suggesting transverse-momentum windows around 3 GeV as promising 
to track QGP radiation. 
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APPENDIX A: PARAMETERISATIONS 



The photon emission rates have been calculated from the Lagrangian describe in Sec.niand by the VMD interaction 

£e„x = -CmlA^pl (Al) 

where is the photon field and C is a constant adjusted by the experimental decay p° — > e+e~, which gives C=0.059. 
In order to respect the Ward Identity in a direct way, we multiply each Feynman amplitude by the square of the 
averaged space- like form factor of Eq. (|10|l . Time-like form factors have been defined to be normalised to one for on- 
shell decays. We quote below parametrisations which include the axial meson oi as exchange particle for non-strange 
initial states. In the following, the photon energy {E) and the temperature (T) are both in GeV. Parameterisations 
for K* K + TT + and K + K^p + j do not appear because their rates have been found to be negligible. 

^dK^^ ^ |^ -(1.46g2°P 0.727) ^ _ o.9957)|^ (fm^^GeV"^) (A2) 



E ^3^^ = F^{E) —exp (^-(9.314r-"-5«4 _ 5.328) (2Ti;)" "«« + (0.3189r"-^^i - 0.8998)-j (A3) 
dR,^^+^+, _ , I f (-35.459rii^6^ 18.827) E\ 



E = ^ (^) T' ^'exp (^- (^^7)1.05 + (2.38 94T" "^«^ - 3.222)- j (A5) 

^ dR^+K^K'+j ^F\E)^exp(-{5.m8T-°''''^*-1.51){2TEf'''-0m-] (A6) 
d'^p T-^ \ T J 



F{E) is the form factor, cf. Sec. Ill Bl 



